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5^ Abstract 

We examine how the qualitative structure of the Doppler peaks in the angular power 
spectrum of the cosmic microwave anisotropy depends on the fundamental nature of the 
perturbations which produced them. The formalism of Hu and Sugiyama is extended to treat 
models with cosmic defects. We discuss how perturbations can be "active" or "passive" and 
"incoherent" or "coherent" , and show how causality and scale invariance play rather different 
roles in these various cases. We find that the existence of secondary Doppler peaks and the 
rough placing of the primary peak unambiguously reflect these basic properties. 
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The cosmic microwave background (CMB) promises to become one of the most successful 
bridges between theory and experiment in cosmology. As the body of experimental data continues 
to grow[l] theorists are evaluating the impact of this data on the two major paradigms for cosmic 
structure formation: inflation [2], and topological defects [3]. One of the important theoretical tools 
is the formalism of Hu and Sugiyama (HS)[4] which we extend here to accommodate topological 
defects. 

The so-called Doppler peaks, in particular, have attracted great interest. They consist of a 
system of oscillations, known to be present for most inflationary models, in the CMB angular 
power spectrum C' at 100 < / < 1500. The peaks' height and position can be used to fix 
with some accuracy combinations of parameters left free in infiationary models [5]. Progress on 
defect Doppler peak predictions has been slow (see however [6, 7, 8]). It was suggested in [8] 
that, regardless of the remaining quantitative uncertainties, one could expect dramatic qualitative 
differences between defect and infiationary Doppler peaks. More concretely, it was pointed out 
that the non-existence of secondary peaks is a robust feature of some defect theories resulting from 
the different role played by randomness and causality in these theories. In this Letter we elaborate 
on how general this feature is, and pin down its controlling factors. 

The idea is to focus on the basic assumptions of infiationary and defect theories, isolate the 
seminal contrasting properties, and perform a qualitative analysis of the structure of the Doppler 
peaks purely based on these properties. Infiationary fiuctuations were produced at a remote epoch, 
and were driven far outside the bubble radius by infiation. The evolution of these fiuctuations is 
linear (until gravitational collapse becomes non-linear at late times), and we call these fiuctuations 
"passive". Also, because all scales observed today have been in causal contact since the onset of 
infiation, causality does not strongly constrain the fiuctuations which result. In contrast, defect 
fiuctuations are continuously seeded by defect evolution, which is a non-linear process. We there- 
fore say these are "active" perturbations. Also, the constraints imposed by causality on defect 
formation and evolution are much greater than than those placed on infiationary perturbations. 

The notions of scale invariance and causality have different implications in these two types 
of theory. A scale invariant gauge-invariant potential $ (the Newtonian potential on subhorizon 
scales) with dimensions L^l"^ has a power spectrum -?($) = (I'I'kP) cx in passive theories (the 
Harrison-Zeldovich spectrum). In active theories the most general counterpart to the Harrison- 
Zeldovich spectrum is -?($) = if F^ikri). Moreover, active perturbations are constrained by 
causality, in the form of integral constraints [9, 10], such as those written in terms of the gauge 
dependent energy-momentum pseudo-tensor of [11, 12]. In [13] we show that the density subject 
to the integral constraint can be written in the gauge-invariant form U = a"^ pAx + p'' + ?>hv' . Here, 
a is the scale factor, h = a/a, p (Ay) is the total matter density (density contrast) and the scalar 
defect stress-energy tensor is given by 0oo = p' , Qoi = kiv' , and Qij = p'bij +(kikj — k^6ij . 
Then following [14], on superhorizon scales P(l{) oc k^ for active perturbations. The Einstein 
equations for the scalar gauge-invariant potentials, $ and ^ are [15, 16] 



Since isotropy requires 11^ and pH/k^ (where pH/k^ is simply related to the quadrupole of the 
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photon and neutrino fluctuations) to be constant for small k, the Einstein equations imply that 
scaling active perturbations produce scaling gauge-invariant potentials, which must be white- 
noise on large scales. In particular P(^! — $) = F{kri)Tf' , with -F(O) a non-zero constant. For 
most realistic defects x'^F{x) will then have a single peak, located at a value of x close to the 
location of the peak of the defect compensated structure function [17]. We will see that the place 
and thickness of the peak in x'^F{x) are deciding features for the Doppler peaks induced by active 
perturbations. 

Active perturbations may also differ from inflation in the way "chance" comes into the theory. 
Randomness occurs in inflation only when the initial conditions are set up. Time evolution is 
linear and deterministic, and may be found by evolving all variables from an initial value equal 
to the square root of their initial variances. By squaring the result one obtains the variables' 
variances at any time. Formally this results from unequal time correlators of the form 

(<I>(fc, ry)<I>(fc', ry')) = - k')a{<i{k, ry)),7(<I>(fc, ry')), (3) 

with (t(-) = \/ P{ ). In defect models however, randomness may intervene in the time evolution 
as well as the initial conditions. Although deterministic in principle, the defect network evolves 
as a result of a complicated non-linear process. If there is strong non-linearity, a given mode 
will be "driven" by interactions with the other modes in a way which will force all different-time 
correlators to zero on a time scale characterized by the "coherence time" Tc(k,ri). Physically 
this means that one has to perform a new "random" draw after each coherence time in order to 
construct a defect history [8]. The counterpart to (3) for incoherent perturbations is 

??')) = - k')P{<!>{k, T]), T]' -T]) . (4) 

For |ry' — ry| = |Ary| > Tc(k, rj) we have P(<^(k, rj), Arj) = 0. When convolving P(<^(k, rj), Arj) with 
functions which vary slowly at the scale of Tc(k, rj) we may implement an approximation where 

ri)(S>(k', ri')) = 6(k - k')6(ri - ri')P''mk, rj)), (5) 

in which 

P''mk,T])) = j dAT]Pmk,T]),AT]) (6) 

is the time-integrated power spectrum [17]. We shall label as coherent and incoherent the per- 
turbations satisfying (3) and (5) respectively. This feature changes the way the average C' are 
computed, resulting in a striking qualitative difference in the structure of Doppler peaks. We 
expect that equation (3) and (5) will be only rough approximations for some defect cases but still 
allow some intuition to be gained. 

A large class of theories is embraced by combinations of the two concepts just introduced. 
Inflationary perturbations are passive coherent perturbations. Defect perturbations are active 
perturbations more or less incoherent depending on the defect. We will submit this class of 
theories to the computational machinery of Hu and Sugiyama [4], which was initially tailored for 
passive coherent perturbations. We consider the limit where SI = 1 and Slj = 0. Generalization 
is straightforward from HS. Whereas primordial terms are dominant for passive perturbations, 
for active perturbations one may drop all but the convolution terms. The radiation brightness 
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multipoles at epoch rjo (rjo > rf , the epoch of decoupling) are then given by 
e'(fc,ryo) = r dvk{^-'f){D'(ri) + V'(ri)) 



fVo 

In' 



[(*-$)A'](ry*)+ Hdriii! -^)A'(ri) (7) 

J n* 



with A (rj) = (2/ + l)ji(k(rio — rj)) and the projected monopole and dipole contributions given by 

v3 

V'ik,7j) = iil+l)j,+iiV^Ax)-lj,_iiV^Ax)))cosAx, 

where Ax = k(ri* — ri)/\/3. The average C' are computed in HS assuming (3). For incoherent 
perturbations one has instead 



D' + V'^ ^ 



21+1 ^ 



J dkk'Pi^' -<!>)( + J dke 



(2/ + 1)2 



tt'' I dk j dqk'P^^ -^){^^] (8) 
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resulting from squaring and averaging (7) using statistics as in (5). These generalizations, together 
with the active perturbations' gauge-invariant potentials, allow a systematic extension of the HS 
formalism to the defect case. 

We first undertake a preliminary analysis by examining the power spectrum of the radiation 
energy density (the monopole term Go + ^) at last scattering. This should mimic the Doppler 
peaks' qualitative structure. As in [4] we have for adiabatic and isocurvature passive fiuctuations: 

k^P(eo + 'S)(r]*) K k^P (So + 'S)(0)cos'^x*, (9) 
k^P(eo + 'S)(r]*) K k^P(eo + 'S)(0)sm'^ X*, (10) 

in which x* = krj* The peaks of the spectrum are at the scales x'^ = rmr for adiabatic, and 
x"^ = (m — l/2)7r for isocurvature perturbations. These correspond roughly to the angular scale 
Im ~ ^mVo- For coherent active perturbations k^P(Qo + ^)(ry*) is approximately 



—= / dx x~ F~ (x) sin x* 



whereas for incoherent perturbations one has 



x*^Fix*) + ^l^ dxx'^Fix)sm'^ (^x* - 

where _P($ — ^) = rpF^x) (or P'^(<^ — xp) = ri'^F(x)), with x = kr/. These show that the position 
and structure of active perturbations Doppler peaks result from a combination of the issue of 
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Figure 1: Three potential structure functions F(x)x^ and their corresponding spectra k^P(Qo + 
assuming coherence. By increasing Xg one moves from adiabatic peaks (line) to isocurvature 
peaks (dash). For larger Xg the secondary peaks come out significantly distorted (point-dash). 

coherence and the details of the potential structure function F(x). If x^F(x) has a sufficiently 
thin peak at x = Xg = 27rry/^c (where is approximately the coherence length of the defect) 
then the monopole spectrum peaks will be at x'^ = (m — l/2)7r + Xc/\/^ for both coherent and 
incoherent fiuctuations. Then active perturbations just apply a phase shift of value Xc/\/3 — 7r/2 
to an adiabatic type of spectrum. 

For Xc 2.7 (unrealistic because it is very close to the smallest turnover point allowed by 
causality[10]) the monopole peaks are at the adiabatic positions. For all other causal active 
perturbations the peaks are shifted to smaller scales. For Xg 3.4 they are out of phase with the 
adiabatic peaks (as in [6]). For Xg > 5.5 the peaks start only in the adiabatic secondary peaks 
region. For standard values of Slj and h these three cases would place the main "Doppler peak" 
at / 230, 350, and 500, respectively. Therefore the placing of the peaks is not a generic feature 
of active fiuctuations. Active perturbations simply add an extra parameter on which the Doppler 
peaks position is strongly dependent. In general we should expect that for the same SI, Slj, and h, 
active perturbations will take the predicted CDM adiabatic peak position / to / + rioixc/\/^ — 7r/2). 
The secondary peaks' separation is not changed, in a first approximation. This is to be contrasted 
with non-fiat infiationary models where C;(S1 = 1) is taken into Ci^i/2. Thus it should be possible 
to distinguish between low SI CDM and = I high-a;c defects. 

For large Xg the peak in x^F(x) can never be thin. Then each mode is active for several 
expansion times, bringing coherence into play. Qualitative changes come about in the secondary 
peaks, but our conclusions relating to the primary peak still hold. In Fig. 1 we consider coherent 
perturbations with realistic structure functions. One may obtain passive type of peaks at adiabatic 
and isocurvature positions. For Xg > b there are strong distortions. One must however realize that 
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Figure 2: The k^P(Qo + ^)(ry*) spectrum for coherent (line) and incoherent (dashed) active 
perturbations with the first two structure functions used in Fig. 1. One may obtain (softer) 
secondary oscillations at the adiabatic position for incoherent perturbations. As the spectrum 
shifts to the right (larger Xg) the secondary oscillations disappear very quickly. 



effective coherence for large Xg requires actual coherence over several (more than Xg) expansion 
times, perhaps an unreasonable demand. 

The situation for incoherent oscillations is illustrated in Fig. 2. Although there are never 
true zeros in P(Qo + it is possible to obtain significant oscillations if the main peak is 

at the adiabatic position. However these disappear very quickly as the main peak approaches 
the isocurvature position. Hence, large x^ defects can always be expected to produce exotic C' 
spectra, as suggested in [8]. Coherence intervenes in deciding how large Xg must be for this to 
happen, as well as what type of novelty is introduced. 

Besides this qualitative general analysis, the extended HS formalism allows for an approximate 
solution (5-10%) for the C''s of any particular defect model. The calculation errors are in practice 
dominated by uncertainties in the defect stress-energy tensor. We illustrate the procedure with 
the example of cosmic strings. For these we use the incoherent form with P'^{p') = 1/(1 -|- 2(j3x)^) 
(from [17]). We consider the two cases /3 = 1 and /3 = .3 similar to the X and I models in [17]. We 
consider only scalar contributions. We assume that the defect variables are subject to equations of 
state of the form = j(x)p'' , H^ = Tf'^s{x)p' , and = ri^v{x)p' ■ Energy conservation at small 
X requires that 87(0) = (l/2a) — 1 and 7t,(0) = (1 — 2a)/(3a(4a -|- 1)), with a = rjh. We make 
use of a string simulation to determine the large x behavior[13]. We find, with large uncertainties, 
that x^y{x) = s .1 — .3, and x'^^'{x) = u k, A — .55. We interpolate between the x ^ 1 
and X <^ 1 behaviour. We set H = and assume that a'^pAx is subdominant except for the 
compensation. This consists of a white noise large-scale tail in the spectrum of a^/oAy present in 
order to cancel the white-noise tail in + 3hv'' and ensure that P(a'^pAT + p'' + 3hv'') goes like k^. 
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Figure 3: The C' spectrum for I (dash) and X (line) cosmic strings. The top lines use s = A, a = A 
for X strings and s = .15, a = .45 for I strings. The bottom lines both use s = .2, a = .5. We have 
assumed SI = 1, /i = .5, and Slj = 0.05. 

We take the compensation into account by hand, setting a^/oAy + + S/it;^ = 7c (p^ + S/it;^) with 
= I / (I + (a / x)^) . We fix the compensation scale at a = 27r. Using (1) and (2) we finally obtain 
the required cosmic strings potential structure functions to be inserted in the HS formalism as 
modified for incoherent perturbations. The results are plotted in Fig. 3. The Sachs- Wolfe plateau 
exhibits a "running" tilt ranging from n 1.4 before / = 10 to n 1.2 at 30 — 40. There is a single 
Doppler bump located at / 400 — 600 These features are remarkably robust against uncertainties 
in the equations of state. The ratio between the peak and the plateau heights, on the other hand, 
can change by as much as an order of magnitude. 

In general, we find that features of the C' spectrum suggested by the monopole spectrum at 
decoupling are confirmed. Generic defects place the main peak to the right of the CDM adiabatic 
peak. Coherent defects exhibit shifted CDM-type secondary oscillations up to the isocurvature 
positions (which are easily distinguished from the shifts associated with varying SI). From then 
on coherent defects show a peculiar type of secondary oscillations. Incoherent defects erase the 
secondary oscillations if the main peak is placed on or to the right of the isocurvature position. 
Thus the most dramatic effects occur for large Xg defects (such as cosmic strings) where the 
C' spectrum shape at 100 < / < 1500 is radically different according to the active/passive, 
and coherent/incoherent nature of the perturbations. The signature becomes progressively less 
prominent pushed to the lower limit imposed by causality. 
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